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Abstract 

We investigate the cosmological production of gravitational waves in a non¬ 
singular flat cosmology powered by a “running vacuum” energy density 
described by p\ = pa{H), a phenomenological expression potentially linked 
with the renormalization group approach in quantum held theory in curved 
spacetimes. The model can be interpreted as a particular case of the class 
recently discussed by Perico et ah (Phys. Rev. D 88, 063531, 2013) which 
is termed complete in the sense that the cosmic evolution occurs between 
two extreme de Sitter stages (early and late time de Sitter phases). The 
gravitational wave equation is derived and its time-dependent part numer¬ 
ically integrated since the primordial de Sitter stage. The generated spec¬ 
trum of gravitons is also compared with the standard calculations where an 
abrupt transition, from the early de Sitter to the radiation phase, is usually 
assumed. It is found that the stochastic background of gravitons is very 
similar to the one predicted by the cosmic concordance model plus inhation 
except at higher frequencies {u > 100 kHz). This remarkable signature of a 
“running vacuum” cosmology combined with the proposed high frequency 
gravitational wave detectors and measurements of the CMB polarization (B- 
modes) may provide a new window to confront more conventional models 
of inflation. 
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1. Introduction 


In the last decade, many authors have proposed cosmological models 
driven by a “running vacuum” energy density, pvac = A{H)/87rG. The 
leitmotiv of such an idea is related to two basic difficulties of the stan¬ 
dard ACDM model (constant vacuum energy density). Firstly, a dynamical 
A(Fr)-term may solve the so-called cosmological constant problem because 
in this case the vacuum energy density relaxes to its present value and one 
may argue that A is small today because the expanding Universe is too 
old. Secondly, a dynamical-A(hf) term may also solve the so-called cos¬ 
mic coincidence problem, i.e. the fact that the time-varying matter energy 
density and the (constant) vacuum energy density have the same order of 
magnitude nowadays. 

Since long ago, many different phenomenological decay laws for A{H) 
were proposed in the literature (see [1, 2, 3, 4, 5] for the oldest literature and 
[6, 7, 8, 9] for more recent articles). The predictions of the latest models 
have also been confronted with the available observational data and the 
results compared with the predictions of the standard ACDM cosmology 
[ 10 , 11 ]. 

Besides enlarging the standard view of a rigid A-term, there are also 
some attempts to justify theoretically a “running vacuum” cosmology based 
on different methods, among them: the thermal instability of a de Sit¬ 
ter spacetime [12] and the renormalization group (RG) approach in curved 
spacetimes [13]. In the latter case, for instance, the emerging dynamical 
A(/7)-term (beyond the true constant vacuum contribution) depends on an 
expansion power series on the Hubble parameter (^ 77”) where only the 
even powers of H are involved in the RG realization, selected by the general 
covariance of the effective action appearing in the quantum field theoretical 
treatment in curved spacetimes [13, 14, 15]. In the same vein, several 
alternative approaches also try to represent the interacting A(77) models 
through a Lagrangian description as in the F{T) gravity approach [16], 
or in its generalized form F{R,T), as discussed by Harko and coworkers 
[17]. Other attempts involve a mixture of a scalar field interacting with 
the radiation bath [18, 19], as in the so-called warm inflationary model 
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[20, 21, 22, 23], or still based on the quantum mechanics probability of un¬ 
stable states [24], Although interesting and highly promising to understand 
the decaying vacuum problem in the evolving Universe, none of them can at 
present be considered definitive and/or widely accepted by the community. 

In this connection, a large class of nonsingular models, where the vac¬ 
uum energy density evolves phenomenologically as a truncated power-series 
in the Hubble parameter, has been proposed a couple of years ago [25, 26] 
(its dominant term Pvac(-ff) cc H'^,n > 2). This class of models has some 
interesting features, among them: a nonsingular origin for the expanding 
universe (no horizon problem) with a deflationary process also without “exit 
problem”, that is, the model evolves smoothly from the primeval nonsin¬ 
gular de Sitter state to the radiation phase; its late-time cosmic expansion 
history is very close to the concordance model and it also furnishes a smooth 
link between the initial and final de Sitter stages through the radiation and 
matter dominated phases. The temperature behavior and the entropy gen¬ 
eration during the continuous non-adiabatic transition from de Sitter to the 
radiation phase has also been investigated [27, 28] and a comparison with 
the late time observations has also been carried out in detail by Gomez- 
Valent and Sola [29]. 

Furthermore, the recent LIGO detection of gravitational waves (GWs) 
opened a new window of observations for astronomy and cosmological prob¬ 
lems [30]. Its importance goes much far beyond the search for signatures of 
compact binary coalescence (black holes, neutron stars, etc). For instance, 
tensor perturbations describing GWS generated in the inflationary stage 
may affect the pattern of Gosmic Background Radiation (GMB) anisotropies 
through the B-modes polarization [31]. In addition, there are also ongoing 
projects (BICEP3, Keck Array experiments, SPT) and future probes with 
high sensitivity instruments like the QUBIG, especially designed to measure 
the B-modes with high precision [32, 33]. Naturally, as compared with more 
conventional inflationary models, the amplihcation of GWs from a primor¬ 
dial de Sitter stage supported by a decaying-A(i7) model may affect GMB 
polarization in a different way, and, as such, this line of inquire deserves a 
closer scrutiny. 

In this context, we analyze here the production of primordial GWs for a 
class of nonsingular “running vacuum cosmologies”. The phenomenological 
A(i7)-term adopted here is defined by: A(i7) = A;, -|- aH^/Hj, where Af, is 
the constant bare vacuum energy density, H is the Hubble parameter, Hj 
is the primordial inflationary scale, and a a dimensionless free parameter. 
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Therefore, unlike the general class discussed in Refs. [25, 26], the very 
late time behavior of the model analized here is exactly ACDM. The GW 
equation is derived and its time-dependent part numerically integrated since 
the primordial de Sitter stage. 

As we shall see, for higher frequencies (z/ > 100 kHz) the predicted 
spectrum departs from the standard inflationary prediction, and, as such, 
these models are distinguishable if high frequency GW detectors become 
operative and reach the expected sensitivity in the near future. In gen¬ 
eral grounds, this signal reinforces the possibility of a new observational 
approach to inflationary physics, and provides additional motivation in the 
search for stochastic background of GWs at high frequencies [34, 35, 36]). It 
is also argued that similar results remain valid for generic decaying vacuum 
cosmologies with an initial de Sitter stage supported by a power-law AT" 
with n > 2. 

2. The model: basic equations 
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Figure 1: (a) Evolution of the quantity R = a'/a during the transition from primordial 
de-Sitter stage to the standard radiation phase (w = 1/3). The blue line curve shows 
the evolution of R in the standard approach (inflation -|- radiation) assuming an abrupt 
transition (AT). The red lines display the evolution of R in the smooth transition case 
(ST) for different values of the initial Hubble parameter Hi defining the early de Sitter 
phase, (b) The same as Fig. la but now for the potential, V = a"/a. Again, the blue 
curve describes the “Grishchuk potential” usually assumed in the AT treatment [40, 41]. 
Note that the maximum height of the potential for the ST case is strongly dependent 
on the values of Hj. To better visualize the plots we have defined a suitable timescale 
Tirf) (see text). The time Tmax' for which the maximum value of % and V are attained 
assuming AT is always delayed in comparison to the ST case discussed here. 
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Let us now consider that the Universe is described by a flat Friedmann- 
Lemaitre-Robertson-Walker (FLRW) geometry. In the co-moving coordi¬ 
nate system, the background line element reads (c = 1): 

ds^ = dt^ — a^{t)df, ( 1 ) 

where a(t) is the scale factor. 

The Einstein equations in the above background read 

87rGp + A{t) = 3H\ (2) 

87rGp-A{t) = -2H-3H\ (3) 

where a dot means time derivative and iL = a/a is the Hubble parameter. 
To solve the above set of equations one needs the functional form of A(t), 
or equivalently, A{H), as well as an equation of state (EoS). 

As remarked before, it will be assumed here that the A(iL)-term is given 
by a particular case of the general class studied by Perico et al. [25, 26] (for 
closely related works see also [3, 4]) 

87rGp,UH) = A{H) = A, + 3a^, (4) 

where Ai, is the bare cosmological constant. By assuming that the material 
medium obeys the EoS, p = up, where cu is a different constant for each 
era, one may show that the scale factor and the Hubble parameter obey the 
following equations (from now on, without loss of generality we consider 
a = 1 ): 


2H -|- 3(1 -|- u^H 



H 


(1 - 1 - u)Af, — 0 . 


( 5 ) 


The standard cosmic concordance model equations are recovered by taking 
the limits Hj » H and a; = 0 (ACDM). At early times, the last term 
proportional to the bare A;, can safely be neglected. In this case, for {u = 
1/3) the solution of the above equation takes the following form: 


1 + Ca^ ’ 

where G is an integration constant. We see that the transition from an early 
de Sitter (Ua^ << 1, iL ~ Hj) to the radiation phase (Ca^ >> 1, a oc 
is analytically described. 
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Now, in terms of the conformal time [dt = a{ri)dri], the line element (1) 
becomes: 

ds"^ = a‘^{r])[dr]^ — dijdx^dx^], (7) 

while the equation of motion (5) takes the form 


2 


H' 


a 


+ 3(1 + uj^H 



H 


(1 + Ci;)Ab — 0, 


( 8 ) 


where primes denote derivative with respect to r]. By hxing the EoS pa¬ 
rameters, the integration constants for each era are obtained by imposing 
the continuity conditions for 0 ( 77 ) and a'( 77 ) at the transition times between 
two subsequent eras. 

By assuming that the vacuum decay mainly into ultra-relativistic par¬ 
ticles {u = 1/3) when the bare term of the decaying vacuum is negligible 
(Afe << H^/Hj), it is easy to see that ( 8 ) boils down to a" — = 0 . 

A direct integration of this equation yields the following solution: 


0 ( 77 ) 


1 


77 + C*2 -|- 


iv + C2y + 


4Ci 


(9) 


which has two limiting cases: at very early times it is a de Sitter solution, 
a oc | 77 |“^, whereas at late time {t >> Hj^) , it enters in the radiation 
phase, a oc 77 [25, 28]. The reduced Hubble parameter, 'Hijf) = a'/a, for this 
stage 


7^Inf-R(?7) = ~ = 


(77 -|- C2Y + 


4Ci 

Ih 


- 1/2 


( 10 ) 


assumes its maximum value for 77 = —C 2 , so 'H(—C 2 ) = = "Hmax = 

Therefore, we can rewrite the integrating constants in the form 
Cl = and C 2 = —T/max, and, naturally, ( 10 ) can be rewritten in 


terms of the new pair of constants {rj^ 
in the next section. 


c'Hmax), a form that will be useful 


3. Cosmological tensor perturbations 

The classical tensor metric perturbation for the conformal FLRW flat 
metric (7) can be written as [39]: 

= 0^(77) [d77^ — (5ij + hij)dx^dx^], 
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( 11 ) 












where the perturbation hij is small, \hij\ <C 1, is transverse-traceless and 
satisfy the gauge constraints: = 0, h- = 0, V^hij = 0. 

The hrst order evolution equation of hi is given by [39, 37, 38] 

hf' + 2-hf - = leTTGa^^f] = 0, (12) 

where in the last equality we used that = 0 because the total EMT 

w/ (i J 

(matter plus vacuum) has the perfect fluid isotropic form [39]. The general 
solution of the above equation can be Fourier expanded in the standard way 




y/lOnG 

(27r)3/2 



X 




Cn + K{ri)e- 



(13) 


' T* 

where hn{T]) are the mode functions, n is the comoving wave vector, Cn and 

Cn are complex numbers, and ejj(n) is the symmetric, transverse-traceless 
polarization tensor. Now, by recalling that the comoving wave number is 
given by 


n=\n\ = ‘^^^^ = ka{r]), (14) 

A 

and inserting the solution (13) into (12), it is readily seen that the temporal 
part yields a differential equation valid for both polarizations (henceforth 
we drop the index r) 

hnir])" + 2'H{r])hn{'ny + n^hniv) = 0- (15) 

where we have conveniently used the reduced Hubble parameter, "H, whose 
continuity also guarantees the continuity of the hrst derivative of the scale 
factor. By using the auxiliary function, p = hn{vi)a{r}), the above equation 
assumes the form hrst derived by Grishchuck [40, 41]: 

)h = 0, (16) 

This is the basic equation which allow us to obtain the associated physical 
quantities like the wave amplitude, energy density and power spectrum. 




n 


a 

a 
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The scale of the “Grishchuck potential”, V{ri) = a"/a, as compared 
with the wave-number determines the behavior of the limiting solutions 
for /i(? 7 ). The solutions are oscillatory, /i oc when ^ \V\ holds. 

The high-frequency modes in this case are diluted by the cosmic expansion 
since h = In the opposite regime one finds p oc a so that the 

amplitude remains constant. In this case, the damping of the waves due the 
universe expansion is avoided, a phenomenon usually referred to as adiabatic 
amplification [40, 41]. 

In Figures la and lb we show the behavior of the quantities 'H{ri) and 
V{ri) during the transition from the early de Sitter stage to the radiation 
phase. For the sake of clarity, both quantities were expressed in terms of 
a convenient variable, r = 'Hmax(hmax — h)) where r^^ax is the time when 
l-L = 'Hmax is reached. As shown below Eq. (18), the value of "Hmax depends 
on the scale Hi. In the variable r, the maximum of H for the decaying 
vacuum model occurs at Tmax = 0 while for the abrupt case the maximum is 
delayed (r^iax > T ma x). The interesting point here is that such a transition 
in our model is smooth and can analytically be followed in terms of the 
scale factor and the quantity H. 


4 . Gravitational wave solutions 


Let us now determine the generated spectrum of the GWs. During the 
transition from inflation to radiation (Inf-R), the GW equation (16) assumes 
the form 

Minf-R + ~ kinf-R(h)) /^Inf-R = 0, (17) 

where Rinf-R(^7) = a"/a = From section 2 one finds that the 

quantity 'Hinf-R for w = 1/3 reads: 


77lnf-R(?7) 


1 


(18) 


where we have at rjmax the maximum value "Hmax — Ffo \/VeqHi /8 and r^eq 
is the conformal time of equality of the energy densities of radiation and 
matter. Since the scenario starts as a de Sitter spacetime (see discussion 
below Eq. (6)), it is possible to obtain an analytical solution for pinf at 
that time. In addition, by imposing the adiabatic vacuum hm„_j.oo yUinf oc 
Constraint [27], we find: 


^0 



Q-ini'n-rima.x) 


Aiinf(n, r]) 


rs> 


n{V - 


(19) 







loglflrms^. 'loW 



Figure 2: (a) Present day root-mean-square amplitude of the GWs as a function of the 
physical frequency for Hj = 10^® s“^. The amplitudes are displayed both for the decaying 
vacuum cosmology (red-curves) as well as for the ACDM model (blue-curves). In the 
former case the model evolves smoothly from inflation to radiation while in the latter 
an abrupt transition is assumed (see also Fig.l). The superposition suggests that the 
predictions are quite similar in the domain of low frequencies. The infographic display 
the difference for high frequencies (z/ > 100 kHz), (b) Density parameter for GWs in the 
domain of high frequencies. For the ACDM model we have fixed the Hubble parameter 
to be Hinf = 10^® s“^ while for the decaying vacuum cosmology the arbitrary scale Hi 
assumed three possible values as indicated in the figure (see text). 


where Aq is a real constant which specihes the initial amplitude of the 
GWs. The above expression for /iinf(n., rj), as well as its derivative rj), 

provide the initial conditions (taken at the initial unstable de Sitter era) for 
the numerical integration of Eq. (17). Consequently, the evolution of the 
GWs can be traced until a time rjr from which the scale factor and waveforms 
behave as in the usual ACDM scenario. At the radiation era, it is possible 
to hnd an analytical solution for whose initial conditions (evaluated at 
?7r) are given by the previous numerical integration. This solution is valid 
until the conformal time ? 7 eq, for which matter {uj = 0) starts to dominate, 
and the scale factor evolves exactly as in the standard matter dominated 
cosmological scenario. At this stage, the analytical solution was obtained by 
considering the initial conditions (at r^eq) given at the end of the precedent 
radiation era. Finally, we have evaluated //(r/o) at the present time r]Q. In 
addition, from the standard dehnition, we also hnd the resulting present 
power spectrum of the relic GWs in this model: 

TT 
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'P(n,ria) 


( 20 ) 



















































In Figure 2a we show the present day root-mean-square (rnis) amplitude 
of the GWs as a function of the physical frequency v for the decaying vacuum 
model, which is related to the power spectrum via hrms{^, Vo) = Vo)- 

For comparison, we also show the rms amplitude for the abrupt three phase 
transition model (de Sitter, radiation and matter) with Hj = 10^® s“^. For 
simplicity, we have not included a late time A;, dominated epoch since the 
same effect for both models is obtained, namely, a little smaller value of 
hrms for the complete GW spectrum [42], 

Notice the remarkable superposition of both spectra for almost the entire 
spectral range. The models are distinguishable only at very high frequencies, 
which are displayed in detail in Figure 2b, where we show the spectrum of 
the density parameter flowho)- In this hgure, we have hxed the value of 
Hj for the abrupt transition model and considered some possible values of 
Hj for the decaying vacuum cosmology. 

At high frequencies the two models predict distinct spectra for a given 
value of H[. As shown in Figure 1, such an effect can be understood in 
terms of the behavior of 'H(h) and V{ri). Since the decaying vacuum model 
evolves smoothly from a de Sitter towards a radiation era, the shape of 'H(h) 
for the transition and the consequent lower value of T^max (for the same Hi) 
result in a lower high frequency GW production. 

On the other hand, there is no adiabatic amplihcation for frequencies v > 
^max so that we have introduced a cut-off at r'max = "Hmax in these graphs. 
For a given Hj value, the cut-off frequency for the abrupt transition model 
is twice the cut-off frequency for the decaying vacuum model. However, 
since we have freedom to choose Hj for the A(H)-model, different cut-off 
frequencies are allowed (see Fig. 2b). In principle, even for different values 
of Hj, the shape of the curves predicted by the AT and ST models are quite 
different thereby suggesting a possible test for the underlying inflationary 
mechanisms. 

5. Final comments 

We have investigated the production of GWs in the context of a flat 
nonsingular decaying vacuum cosmology. The dynamical A(H)-term was 
phenomenologically modeled as A{H) = Af, -|- H^/Hj. This kind of model 
undergone a smooth transition from an early inflation to the standard radi¬ 
ation phase which can analytically be described [3, 4, 25, 26]. Interestingly, 
the model explain the present day entropy content of the Universe stored 
in the relic radiation [27, 28] without necessity of a highly non-adiabatic 
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reheating stage. Such scenario is also free of horizon and “graceful” exit 
problems, and is also in agreement with the observations at low and in¬ 
termediate redshifts since it evolves to the standard A^CDM cosmology. 
In comparison with the standard approach (big-bang, adiabatic inflation, 
reheating, plus radiation and subsequente eras), the present study also pro¬ 
vides a simple way to understand how a different scenario can modify the 
predictions concerning the generated spectrum of GWs. 

Solutions for the GW equation were numerically obtained. The result 
furnishes a simple and dehnite example that details of the transition from 
an early de Sitter to the radiation phase plays an important role in the 
generation of the GW spectrum. The identihed lower GW production for 
high frequencies (z/ > 100 kHz) is a remarkable signature of the A decay 
model (see Figs, lb and 2b). Although far exceeding the frequencies (and 
sensitivity) of the existing GW detectors in ground-based experiments and 
space-based successors like LISA, new kinds of interferometers and detec¬ 
tors operating in a frequency range high enough to test primordial GW 
have been discussed in the recent literature [43] (see also [44, 45]). The high 
frequency behavior shows that the model is not only distinguishable from 
abrupt inflationary scenario but can also provide a crucial test for the un¬ 
derlying mechanism as predicted in some multi-held inhationary models [35]. 
Note also that current GW interferometers can not detect the putative GW 
background studied here, since the amplitude of this background is orders of 
magnitude lower then the interferometer sensitivities. We refer the reader 
to the paper by Moore, Gole & Berry [46], where the sensitivity curves of 
several detectors are displayed. Finally, it should be reinforced that the 
model discussed here is a starting point for the investigation of more com¬ 
plex and rich decaying vacuum cosmologies, like the general class proposed 
in Ref. [25, 26]. Such models with A{H) oc H'^,n > 2 deserves a closer 
scrutiny since they furnish a complete cosmological history. In principle, al¬ 
though the details might be dependent on the power n, similar results to the 
ones derived here should be expected because the entire class also evolves 
through a smooth transition from de Sitter to the radiation. A detailed 
analysis involving the generation of GWs in this more general framework 
plus the implications in the observed pattern of GMB anisotropies through 
the B-modes polarization will be discussed in a forthcoming communication. 
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